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FOREWORD 


This  report  was  prepared  by  MITHRAS,  Inc.  ,  of  Cambridge, 
Massachusetts,  for  the  Upper  Atmosphere  Physics  Laboratory,  Iono- 
sphe  ric  Perturbations  Branch,  Air  Force  Cambridge  Research  Labora¬ 
tories,  L.  G.  Hanscom  Field,  Bedford,  Massachusetts  under  Contract 
AF 1 9(628)  -42 18. 

The  investigations  whose  results  are  reported  here  were  con¬ 
ducted  during  the  period  1  July  196  5  to  3  May  1966.  The  final  report 
was  written  by  Mr.  P.  O.  Jarvinen.  The  addendum  to  this  report,  which 
deals  with  an  ionization  model  for  high  altitude  rocket  plumes  is  classi¬ 
fied  SECRET,  Request  for  copies  of  the  addendum  should  be  submitted 
to 

Mr.  William  F.  Ring  CRUR 

Air  Force  Cambridge  Research  Laboratories 

Laurence  G.  Hanscom  Field 

Bedford,  Massachusetts 

This  report  concludes  the  work  on  Contract  AF  19(628) -4218. 
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ABSTRACT 


The  gasdynamic  structure  of  high  altitude  rocket  plumes  is  in¬ 
vestigated.  A  3imple  analytical  model  for  the  expansion  of  a  gas  into  a 
vacuum  is  constructed  and  is  shown  to  represent  the  expansion  of  ex¬ 
haust  gases  from  rocket  motors  into  the  region  bounded  by  the  inner 
shock  wave  and  Mach  disc.  A  comparison  of  the  exhaust  expansion 
model  with  solutions  obtained  with  the  method  of  characteristics  is  made 
and  shows  good  agreement.  Contours  of  c  onstant  density ,  electron  den¬ 
sity,  and  collision  frequency  in  the  exhaust  plume  of  a  typical  rocket 
engine  are  determined.  It  is  noted  that  the  theory  of  Alden  et  al,  which 
determined  the  geometric  features  of  the  plume  surface  (i.  e.  contact 
surface)  maybe  used  to  define  the  volume  in  which  exhaust  flow  proper¬ 
ties  are  predicted  by  the  exhaust  expansion  model.  The  plume  scaling 
parameter  along  a  typical  trajectory  is  noted. 

The  blast  wave  theory  of  Hill  et  al  is  extended  to  allow  the  cal¬ 
culation  of  the  size  and  shape  of  rocket  plumes  of  solid  propellant  mis¬ 
siles.  A  method  is  derived  which  determines  the  plume  drag  of  solid 
propellant  motors.  The  method  accounts  for  the  fact  that  solid  propel¬ 
lant  rocket  plume  flow  is  a  two  phase  flow  of  gas  and  solids  and  that  the 
solids  present  in  the  exhaust  take  part  in  determining  the  mass  flow  but 
not  the  pressure  or  the  expansion  of  the  exhaust  gases  upon  exit  from 
the  rocket  nozzle.  The  plume  expansion  and  resulting  plume  size  is 
shown  to  depend  on  the  thrust  and  the  engine  exit  plane  conditions  of  the 
gas  phase  of  the  two  phase  flow. 

A  method  is  described  which  allows  the  average  prftperties  in 
the  air -exhaust  gas  layer  on  the  periphery  of  the  olume  to  be  evaluated. 
The  variation  of  mass  flow,  velocity,  pressure,  temperature  and  over¬ 
all  thickness  along  the  layer  is  ascertained.  The  analysis  is  extended 
to  give  the  average  properties  of  the  air  layer  between  the  dividing 
streamline  and  the  outer  shock  and  the  exhaust  gas  layer  on  the  inside  . 
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The  theory  of  Whitham  (1952)  which  describes  the  far-field  flow  about 
supersonic  axisymmetric  bodies,  is  used  to  evaluate  the  plume  drag 
and  is  found  to  be  in  agreement  with  the  blast  wave  predictions  of  Hill 
et  al. 

In  the  addendum  to  this  report,  a  high  altitude  plume  ionization 
model  is  developed  using  the  theoretical  models  for  exhaust  flow,  two 
phase  rocket  motor  flow  and  properties  of  the  air-exhaust  gas  layer. 
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1.  INTRODUCTION 


The  aerodynamic  features  in  the  nose  region  of  a  high  altitude 
plume  are  shown  in  Figure  1.  The  rocket  motor  exit  plane  oressure 
is  many  orders  of  magnitude  larger  than  the  local  ambient  pressure. 

The  exhaust  gases  expand  rapidly  after  exit  from  the  rocket  motor, 
reaching  a  maximum  angle  of  expansion  0^  which  is  determined  by  the 
rocket  motor  exit  Mach  number  and  the  ratio  of  specific  heats,  y,  of  the 
exhaust  gases.  The  exhaust  gas  plume  appears  as  a  blunt  body  to  the 
free  stream  flow  and  a  bow  shock  wave  is  formed  about  the  plume. 

The  free  stream  flow  which  passes  through  the  bow  shock,  flows  outward 
in  the  regioi  between  the  bow  shock  and  the  contact  surface.  The  inter¬ 
nal  shock  wave  is  formed  in  a  not  too  obvious  way  by  the  coalescence  of 
compression  waves  from  the  jet  edge.  These  waves  are  simply  reflec¬ 
tions  of  expansion  waves  formed  at  the  exit  of  the  underexpanded  nozzle. 
The  expanding  exhaust  gases  are  deflected  into  the  downstream  direction 
by  the  free  stream  flow  and  are  confined  to  the  region  to  the  right  of  the 
contact  surface.  The  exhaust  gases  which  pass  through  the  inner  shock 
flow  outward  in  the  region  between  the  contact  surface  and  the  inner 
shock  wave  (i.  e.  ,  streamtube  A3).  Exhaust  gases  expand  into  the 
region  bounded  by  the  inner  shock  wave  and  Mach  disc  as  if  into  a 
vacuum.  The  launch  vehicle  is  immersed  in  a  layer  between  the  bow 
and  inner  shocks  and  plays  a  negligible  part  in  determining  the  flow 
characteristics.  The  flow  is  hypersonic  and  continuum  theory  holds. 

Analytical  techniques  have  previously  been  developed  which 
describe  the  inviscid  gasdynamic  structure  of  high  rocket  plumes.  Hill 
et  al  (1963)  developed  a  simple  theory  to  describe  the  size  and  shape  of 
the  plumes  of  expanding  gases  behind  a  missile  in  powered  flight  at 
altitudes  above  100  km  using  the  blasv-wave  analogy.  Expressions  were 
derived  for  the  rocket  p'ume  nose  radius,  maximum  plume  diameter 
and  distance  to  the  Mach  disc.  A  comparison  of  the  theoretical  predic¬ 
tions  with  optical  observations  of  Atlas  and  Titan  plumes  showed  ex¬ 
cellent  agreement  in  the  nose  radii  and  maximum  radii  of  the  plumes. 
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This  theory  relates  the  plume  parameters,  such  as  nose  radius,  to  the 
plume  drag.  The  relationship  between  plume  drag  and  engine  thrust  is 
determined  by  applying  a  momentum  balance  to  the  rocket  plume.  It  is 
shown  that  the  geometric  features  of  the  plume  can  be  described  if  the 
engine  thrust  and  thrust  coefficient,  vehicle  altitude  and  forward  speed 
are  known. 

Alden  et  al,  (1964)  developed  a  simple  inviscid  analytical 
theory  which  defined  the  geometrical  features  of  the  plume  surface 
(i.  e.  ,  contact  surface)  by  applying  a  force-momentum  balance  to  the 
plume  and  determined  a  one -parameter  family  of  functions  which  described 
the  jet  flow  from  rocket  motors.  In  addition,  the  validity  of  using  con¬ 
tinuum  fluid  mechanics  to  describe  high  altitude  plum.;  phenomenon 
was  examined. 

In  this  report  some  additional  theoretical  models  for  the  inviscid 
gasdynamic  structure  of  high  altitude  plumes  are  discussed.  An  improved 
analytic  method  for  calculating  the  expansion  of  gas  from  the  exit  of  a 
rocket  nozzle  into  the  region  bounded  by  the  inner  shock  wave  and  Mach 
disc  is  constructed.  The  method  is  especially  useful  for  high  altitude 
calculations  where  the  plume  extends  for  thousands  of  nozzle  exit  radii 
and  where  the  method  of  characteristics,  while  capable  in  principle  of 
describing  the  exhaust  gas  flow,  has  practical  limitations  such  as  its 
relatively  high  cost  and  the  difficulty  of  calculating  more  than  a  few 
hundred  radii  from  the  nozzle  exit.  A  comparison  of  the  exhaust  ex¬ 
pansion  model  with  solutions  obtained  .with  the  method  of  characteristics 
is  made  and  shows  good  agreement.  Contours  of  constant  density, 
electron  density  and  collision  frequency  in  the  exhaust  plume  of  a  typical 
rocket  engine  are  determined  using  this  exhaust  model. 

The  blast  wave  theory  of  Hill,  et  al,  is  extended  to  allow  the 
calculation  of  the  size  and  shape  of  rocket  plumes  of  solid  propellant 
missiles.  A  method  is  derived  which  determines  the  plume  drag  of 
solid  propellant  motors.  The  method  accounts  for  the  fact  ths.t  solid 
propellant  rocket  plume  flow  is  a  two  phase  flow  of  gas  and  solids  and 
that  the  solids  present  in  the  exhaust  take  part  in  determining  the  mass 
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flow  but  not  the  pressure  or  the  expansion  o£  the  exhaust  gases  upon 
exit  from  the  rocket  nozzle.  The  plume  expansion  and  resulting  plume 
size  is  shown  to  depend  on  the  thrust  and  the  engine  exit  plane  conditions 
of  the  gas  phase  of  the  two  phase  flow. 

A  method  is  described  which  allows  the  average  properties  in 
the  air-exhaust  gas  layer  or  the  periphery  of  the  plume  to  be  evaluated. 
Thf  variation  of  mass  flow,  velocity,  pressure,  temperature  and  over¬ 
all  thickness  along  the  layer  is  ascertained.  The  analysis  is  extended 
1.0  give  the  average  properties  of  the  air  layer  between  tne  dividing 
streamline  and  the  outer  shock  ar,d  the  exhaust  layer  on  the  inside. 
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2.  ANALYTIC  DESCRIPTION  OF  PLUME  FLOWS 


The  expansion  of  a  gas  from  the  exit  of  a  rocket  nozzle  into  a 
vacuum  may  be  calculated  in  a  straightforward  manner  by  the  numerical 
method  of  characteristics.  Although  most  of  the  published  calculations 
consider  the  thermodynamic  behavior  of  the  gas  to  be  perfect,  this 
restriction  is  not  necessary. 

In  most  rocket  plume  problems  of  interest  the  expansion  is  into 
a  finite  back  pressure  rather  than  into  a  vacuum.  A  shock  wave  then 
forms  in  the  jet  near  its  outer  boundary.  The  vacuum  solution  is  useful, 
however,  in  that  it  still  describes  the  flow  inside  the  region  boui.ded  by 
this  shock.  It  may  also  be  used  to  calculate  the  shock  envelope  in  the 
manner  described  by  Alden  et  al. 

Now  the  method  of  characteristics,  while  capable  in  principle 
of  describing  the  flow  from  any  rocket  motor,  has  two  practical  limita¬ 
tions.  One  is  its  relatively  high  cost  and  the  other  the  difficulty  of 
carrying  out  the  calculations  more  than  a  few  hundred  radii  from  the 
nozzle  exit.  Thus  a  simple  analytic  method  may  be  very  useful,  espe¬ 
cially  for  high-altitude  calculations  where  the  plume  extends  for  thousands 
of  nozzle  exit  radii.  Such  a  method,  drawing  on  information  available 
from  exact  solutions,  is  given  below. 

2.  1  General  Character  of  Plume  Flows 

A  schematic  illustration  of  the  flow  exhausting  from  a  nozzle  into 
a  vacuum  is  given  in  Figure  2.  The  expansion  is  bounded  by  the  limiting 
streamline  at  the  angle  ®max»  which  represents  a  Prandtl-Meyer  flow  at 
the  nozzle  lip.  Along  the  other  stream-lines  the  flow  turns  more  grad¬ 
ually  but  "exact*  calculations  show  that  as  the  local  velocity  approacnes 
the  maximum  possible  velocity,  V  ,  all  streamlines  become  essentially 
straight.  This  behavior  is  in  fact  characteristic  of  any  irrotational  flow 
pattern  in  which  the  velocity  variation  from  one  streamline  to  the  next 
tends  to  zero. 
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With  the  velocity  fixed  in  the  far  field  of  the  flow  pattern,  con¬ 
servation  of  mass  along  any  streamtube  requires  an  inverse  square  law 
variation  of  density  with  radial  distance  from  the  exit.  It  is  known, 
m  reover,  that  the  mass  flux  is  a  maximum  along  the  axis  and  falls  off 
very  rapidly  with  angle  of  inclination  from  the  axis.  A  complete  simple 
description  of  the  flow  can  be  given  by  specifying  this  angular  distribu¬ 
tion  of  mass  flux. 


2.  2  Angular  Distribution  of  the  Mass  Efflux 

2.  2.  J  One  Parameter  Family  of  Distribution 

In  our  approximate  description  of  the  flow  far  from  the 
nozzle  the  streamlines  are  straight  and  appear  to  originate  from  a 
common  point  close  to  the  center  of  the  nozzle  exit.  Along  each  stream¬ 
line  we  have 


pVr^  =  constant  =  dm(0) 


(1) 


where 


p  is  the  density 

V  is  the  velocity 

r  is  the  radial  coordinate 

The  quantity  dm  ((?/  represents  the  mass  flux  per  unit  solid  angle  in  the 
direction  0. 

It  is  convenient  to  normalize  dm  (0)  with  respect  to  its  value  on 
the  nozzle  axis  and  to  write 


f(0)  = 

dm  (0) 


(pVrZ)e 

<pVr2)0 


(2) 
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We  then  postulate  that  a  one-parameter  family  of  functions  £(Q)  is  suffi¬ 
cient  to  describe  adequately  the  flow  from  all  nozzles. 

A  suitable  form  of  the  function  i(6)  may  be  obtained  from  some 
available  "exact"  solutions  by  the  method  of  characteristics.  Alden 
et  al  (1964)  has  suggested  the  form 


m  = 


cos  0  -  cos  9 
1  -  cos  9 


max 


max 


m 


where  9  is  the  limiting  flow  inclination  for  expansion  into  a  vacuum. 

TaISLX 

A  somewhat  better  fit  to  a  larger  sample  of  exact  calculations  has  since 
been  obtained  with  the  function 


f (9)  -  exp 


-\'2(1  -  cos  9)Z 


(3) 


A  typical  sample  of  "exact"  results  plotted  in  a  coordinate  system 
appropriate  to  this  function  is  shown  in  Figure  3.  Because  data  on  the 
axis  was  not  available  in  these  cases  we  have  plotted  the  function 

(pVr2)0 

,  ,  ,2  m  =  F(«)  (4) 

par 

The  value  of  the  constant  multiplier  of  f (6)  varies  from  one  nozzle  to  the 
next. 

Unlike  Alden1  s  (1964)  function,  this  new  one  does  not  satisfy  the 
condition 


f  (9 


max 


=  0 


However,  as  shown  in  Figure  3,  *(®inax)  *8  certainly  extremely  small 

and  the  contribution  of  the  fictitious  mass  flux  beyond  9  to  calculations 

7  max 

of  total  mass  and  momentum  fluxes  is  safely  neglected. 
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2.  2.  2  Determination  of  the  Value  of  the  Parameter  X, 

The  parameter  X.  is  determined  by  requiring  that  the 
mass  and  momentum  fluxes  in  the  flow  equal  those  at  the  nozzle  exit. 

These  fluxes  may  be  computed  by  integrating  over  the  solid  angle  con¬ 
taining  the  flow.  We  shall  regard  this  integration  as  being  performed 
over  a  surface  and  in  particular  we  shall  use  constant-property  surfaces 
on  which  the  values  of  the  density,  velocity,  etc.  ,  are  constant. 

The  total  mass  flux  in  the  flow  may  be  written 


°  o 

m  =  2 


■c 


max  ,.2  Q  a 
pvr  sin  0  d  0 


=  Zn  { pVr  ), 


s; 


max 


f (6)  sin  0  d  0 


*  *  *2 

=  Trp  a  r 


(5) 


where  the  superscript  *  refers  to  conditions  at  the  nozzle  throat. 

The  momentum  balance  between  the  constant-property  surface 
and  the  nozzle  exit  plane  may  be  written 


I 


7  .0 

F  =  2t r  (pVr  )QV  |  maX  f(0)  cos  0  sin  0  d  0  +  p  A 


ex 


*CF  Pcr 


*Z 


(6) 


where 


F  is  vacuum  thrust  of  the  nozzle 

A  is  the  nozzle  exit  area 
ex 

V  is  the  velocity  on  the  integration  surface 
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*  * 


p  is  the  pressure  on  the  integration  surface 
Pc  is  the  combustion  chamber  pressure 
Cp  is  the  vacuum  thrust  coefficient 

Now  in  the  far  field  p/pc  is  very  small  and  the  contribution  of  the 
pressure  term  may  be  neglected.  Dividing  the  first  expression  into  the 
second  we  obtain 


s: 


max 


f(0)  cos  9  sin  9  d  9 


s: 


max 


f(9)  sin  9  d  9 


CF  Pc 
p*  a*  V 


(7) 


Now  it  may  be  shown  that: 


P*  a*  V 

11  -  r 

PL  "  CF 


(8) 


max 


where 


is  the  limiting  velocity 


CF  is  the  maximum  vacuum  thrust  coefficient 
max 


Thus 


CF*c 


C.  V 
F  m 


*  Cr 

p  a  V  F 


V  w  C 


Cf\ 


max 


F  \ 
max  \ 


(9) 


The  integrals  are  conveniently  evaluated  in  terms  of  the  param¬ 
eter  X.  by  making  the  substitution 


r|  =  1  -  coe  9 


(10) 
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and,  in  view  of  the  small  value  of  the  integrand  at  ®rnax»  letting  the 
upper  limit  of  t]  -»  co.  Thus 


I 


»  2  2 

e  ^  ( 1  rj)  d  r| 

0 


w  C, 


max 


si 


2  2 
t-*  n  d  ri 


1  ±- 


\*Jn 


(ID 


Choosing  the  sign  to  yield  a  positive  value  for  X  we  finally  obtain 


X 


_1_ 

's/it 


_  1 _ 

1  Cp/w  CF 

max 


(12) 


2,  3  Properties  of  Typical  Flows 


on 


2.  3.  1  Variation  of  X 

The  value  of  the  parameter  X  as  derived  above  depends 


(a)  the  quantity  Cp/Cp  which  is  fixed  by  the  charac- 

£  max 

teristics  of  the  propellant  and  rocket  nozzle 

(b)  the  velocity  ratio  w  which  varies  from  one  constant- 
property  contour  to  the  next 

For  the  asymptotic  flow  at  very  large  distances  from  the  exit 
w  -*•  1  and  the  constant-property  contours  become  self-similar,  fixed 
by  the  parameter  CF/CFmax-  Figure  4  illustrates  the  variation  of  X^, 
the  asymptotic  value  of  X,  with  che  nozzle  area  ratio  and  the  value  of 
y  (ratio  of  specific  heats)  of  the  exhaust  gases.  For  any  given  rocket 
motor  the  value  of  X  varies  along  the  plume.  This  is  illustrated  in 
Figure  5  for  a  set  of  nozzles  with  A/A  =  25. 


2,  3.  2  Density  Decay  Along  the  Axil 


to  yield 


The  mass-flux  balance  (5)  above  may  be  manipulated 


(Pr\ 

*  *  2 

P  r 


(1. 


for  the  variation  of  density  along  the  axis.  Normalizing  the  density  with 
respect  to  its  value  in  the  combustion  chamber, 


(Pr2)0 

*  2 

P  r 
Kc 


Y  +  1 


V  -  1 


w_ 
w 


n/tt 


(14) 


For  a  given  value  of  p/pc»  then  the  axial  station,  (r)^  =  x,  i6  given  by 


x  2 

(JV)  = 

r 


1 


Y+  1 


Y  -  1  .  *  p 

1  X.  w  Kc 


\f T 


W 


(15) 


Thus  the  rate  of  density  decay  depends  primarily  on  the  value  of  \,  but 
is  modified  near  the  nozzle  by  the  variation  in  w.  In  the  asymptotic 
flow  with  w  s  1. 


1 


Y+  1 


* 

w 


(16) 


with  a  constant  of  proportionality  dependent  on  X.^  and  y.  Figure  6  shows 

the  variation  of  density  along  the  axis  for  three  values  of  y.  The  value 

of  X.  in  each  case  has  been  chosen  to  correspond  to  an  area  ratio 
oo 

A  /A*  =  2  5.  Note  that  the  higher  densities  persist  further  with  the 
higher  values  of  y. 
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2.  3.  3  Constant  Density  Contours 


Once  we  have  located  the  point  on  the  axis  corresponding 
to  any  given  density  the  rest  of  the  curve  along  which  the  density  has 
this  value  is  given  by  the  equation 


rZ(9) 

r2(0) 


f  (9)  =  exp 


( 1  -  ccs  9) 


(17) 


that 


x 


Tq  cos  9  exp 


(1 


cos  9)  /2 


(18) 


y  =  Tq  sin  9  exp 


-\2  (1  -  cos  9)2/ 2 


(19) 


Figure  7  illustrates  typical  constant  property  contours  for  the  three 
values  of  \  used  in  Figure  6.  As  expected,  the  narrower  plumes  are 
those  in  which  high  densities  persist  further  downstream. 


2.  4  Comparison  with  Published  Numerical  Solutions 

2.  4.  1  Density  Decay  Along  the  Axis 

In  a  discussion  of  a  different  far -field  approximation  for 
plume  flows,  Sibulkin  and  Gallaher  (1963)  have  published  the  results  of 
some  numerical  calculations  by  the  method  of  characteristics.  They 
describe  the  rate  of  density  decay  along  the  axis  of  the  flow  in  terms  of 
the  parameter 

2 

B  =  -2-  (-^r)  (20) 

Pc  d* 
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evaluated  at  x/re  =  100  for  y  =  1.2,  1.  3  and  1,4  and  at  x/rg  =  50  for 
Y  =  1.  67. 

In  terms  of  the  theory  of  this  report, 


In  Figure  8,  this  formula  is  compared  with  the  computations  quoted  by 
Sibulkin  and  Gallaher  (1963).  The  agreement  is  good.  The  approximate 
theory  errs  on  the  low  side  for  the  low«r  values  of  y  but  by  not  more 
than  30  percent  at  most.  This  error  is  much  smaller  than  that  of 
Sibulkin  and  Gallaher' s  (1963)  method. 

2,4.  2  Constant-Property  Contours 

A  typical  set  of  constant-property  contours  calculated 
numerically  by  the  method  of  characteristics  has  been  published  by 
Altshuler  (19  58).  This  set  is  restricted  to  the  region  near  the  nozzle 
exit  where  the  flow  has  not  yet  become  asymptotic  so  that  the  value  of 
X  changes  appreciably  from  one  contour  to  the  next. 

Figure  9  illustrates  the  three  contours  furthest  from  the  nozzle 
together  with  the  predictions  of  the  approximate  method  of  thi6  report. 
The  agreement  is  good  and  improves  with  increasing  distance  as  ex  - 
pected.  At  the  larger  distances  for  which  this  analytic  method  is  in¬ 
tended  it  should  yield  a  flow  pattern  accurate  enough  for  any  engineering 
application. 

2.  5  Contours  of  Constant  Electron  Density  and  Collision 
_ Frequency  In  a  Typical  Plume. 

The  exhaust  flow  expands  as  if  into  a  vacuum  upon  exit  from 
rocket  motors  operating  at  high  altitudes.  The  rapid  expansion  causes 
the  flow  downstream  of  the  exit  plane  to  be  frozen  (i.  e.  ,  no  chemical 
reactions  occur  to  change  the  exhaust  gas  composition  or  electron 
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population).  Downstream  of  the  engine  exit  plane  the  electron  density 
varies  exactly  as  the  gas  density.  Thus  contours  of  constant  gas  density 
are  also  contours  of  constant  electron  density  (plasma  frequency).  The  same 
contours  also  represent  constant  collision  frequency  contours.  A  knowl¬ 
edge  of  the  gas  density,  electron  density  and  collision  frequency  at  the 
nozzle  exit  plane  in  conjunction  with  the  engine  y  is  sufficient  to  define 
these  parameters  throughout  the  exhaust  flow  within  the  inner  shock  and 
Mach  disc.  Electron  density  and  collision  frequency  contours  have  been 
evaluated  for  a  typical  engine  using  the  analytical  model  for  exhaust  flow 
into  a  vacuum  and  assuming  frozen  flow  downstream  of  the  nozzle  exit 
plane.  (Figure  10).  The  initial  electron  density  and  collision  frequency 
at  the  engine  exit  plane  were  normalized  to  1.  0  electron  per  cubic  centi¬ 
meter  and  1.0  collision  per  second  respectively. 


2.  6  Locations  of  Plume  Surface 


Alden  et  al,  developed  a  simple,  analytic,  inviscid  theory  to 
describe  the  surface  of  rocket  plumes.  The  surface  of  the  plume  is 
described  by  the  dividing  streamline  (contact  surface)  separating  the 
outer  air  flow  from  the  inner  flow  of  rocket  exhaust  gases.  The  location 
of  the  surface  is  determined  by  formulating  a  force -momentum  balance 
normal  to  the  plume  surface  in' conjunction  with  a  model  for  the  exhaust 
flow  field  and  the  assumption  that  hypersonic  flow  theory  is  applicable. 

The  inviscid  plume  surface  is  found  to  be  represented  by  a  nondimensional, 
second  order  non-linear  differential  equation  in  terms  of  r  and  0 
(Figure  11),  the  nondimensional  spherical  radius  and  polar  angle  respec¬ 
tively,  which  when  integrated  located  the  surface  of  the  dividing  stream¬ 
line.  The  differential  equation  contains,  in  addition  to  the  coordinate 
variables  (r,  0),  only  two  explicit  parameters  (y,  0j^),  which  are  fixed 
by  the  rocket  motor  exhaust  properties  and  geometry.  -The  effects 
of  all  other  parameters  (i.  e.  ,  combustion  pressure,  motor  size,  forward 
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speed  and  altitude)  are  combined  impiicity  in  r  which  is  defined  as 


i 


qQ  =  free  stream  dynamic  pressure 
Pc  =  engine  chamber  pressure 
y*  =  engine  throat  radius 


For  a  given  motor  type,  a  single  computation  will  give  the  plume 
shape  for  all  (high)  altitudes,  all  motor  sizes,  various  chamber  pres¬ 
sures,  and  all  hypersonic  flight  speeds.  A  typical  plume  contact  surface 

shape,  calculated  using  this  method,  is  shown  in  Figure  12  for  y.  =  1.  25 

J  Ct 

and  maximum  Prandtl-Meyer  expansion  angles  0^  of  120.  5*,  corre¬ 
sponding  to  exit  Mach  number  of  3. 


Dimensional  contact  surface  geometries  (r  as  a  function  of  6) 
may  be  obtained  from  the  non-dimensional  results  through  the  scaling 
factor 


The  variation  of  the  plume  scaling  parameter  along  a  representative 
trajectory  is  noted  in  Figure  13.  The  size  of  the  plume  is  seen  to  change 
by  one  order  of  magnitude  between  90  and  120  kilometers  and  approxi¬ 
mately  another  order  of  magnitude  between  1.20  and  180  kilometers. 

The  exhaust  flow  within  the  inner  shock  wave  (approximately  represented 
by  the  contact  surface)  remains  unmodified  by  extex'nal  flow  conditions. 
Therefore  superposition  of  the  dimensional  contact  surface  geometry 
on  dimensional  electron  density  contours  specifies  the  volume  in  which 
exhaust  flow  properties  as  predicted  by  the  exhaust  expansion  model 
are  unmodified  by  the  external  flow  field. 
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As  noted  in  the  following  section,  the  plume  surface  for  solid 
propellant  rocket  motors  is  properly  located  with  the  theory  of  Alden 
et_al  if  the  gas  phase  exhaust  pioperties  of  the  two  phase  flow  are  used. 
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The  rocket  nozzle  flow  from  solid  propellant  motors  is  a  two- 

phase  flow  of  gas  and  solids  and  parameters  such  as  C  ,  R  C  and  y 

P  »  v 

must  be  determined  differently  than  for  single  phase  gas  flow. 

Simple  calculations  of  two-phase  nozzle  flows  can  be  made  once 
the  effective  mean  C^,  R  and  y  for  the  gas  phase,  the  effective  mean 
Cp  lor  the  solid  phase  and  the  relative  weight  flows  of  the  gas  and  solid 
pha  ses  are  determined. 

For  a  mixture  of  gases,  the  mole  fraction  n^  is  defined  as  the 
ratio  of  the  number  of  moles  present  of  the  i^  component  to  the  total 
number  of  moles,  both  per  unit  volume.  Thus 


n.  n. 
1  l 

n.  =  -  = — 

1  Zn.  n 

o  1 


The  mass  concentration  is 


M .  n. 
l  l 

r  -  _ 2. 

l  Z.M.n. 
i  ii 


where  is  the  molecular  weight  of  the  i*“  component.  The  specific 
heat  0^,  for  a  mixture  is  the  sum  of  the  component  specific  heats 
weighted  according  to  their  respective  mass  fractions  C.,  as  follows: 

C  =  Z.  C.  C  (25) 

p  i  i  p.  ' 
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The  gas  constant  for  the  mixture  is  determined  from  the  universal  gas 
constant  $  and  the  molecular  weight  of  the  mixture  M 

R  =  S/M  (26) 

The  ratio  of  specific  heats,  y»  f°r  the  mixture  of  gases  is  obtained  from 


C  and  R  as  follows 

P 

C  =  C  -  R 

v  p 

(27) 

Y  =  C  /C 

•  p  V 

(28) 

The  effective  values  of  y  «  R  and  C  for  the  gas -solid  rocket  motor  flow 

L 

are  formed  from  the  following  expressions, 

1  +  K  C  /C 

Ps  P 

Y  Y  1  +  Y  K  C  7C~ 

Ps  P 

(29) 

R 

R  1  +  K 

(30) 

C  +  K  C 
_  P  Ps 

Cp  '  1  +  K 

(31) 

where 

y  mass  flow  of  solids 
mass  flow  of  gas 

(32) 

C  =  effective  specific  heat  of  solids 

D 

s 

(33) 

The  nozzle  flow  conditions  for  the  study  engine  are 

obtained 

from  gas  flow  tables  (Wang  et  al  (19  57))  using  the  engine  expansion 
ratio,  Ae/A  ,  and  the  effective  specific  heat  ratio  y. 
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The  fol’owing  ratios  are  obtained; 


VVM=C1 

Me=C2 

P  /P  =  C, 
e'  c  3 

T  /T  =  C, 
e  c  4 

The  maximum  gas  velocity,  VM>  is  determined  from  (34) 


=  Vc 


1 


(34) 

(35) 

(36) 

(37) 

(38) 


where  Ve-  the  exit  velocity,  is  usually  a  specified  engine  characteristic 
or  is  calculated  for  a  properly  expanded  engine  as 


V 

e 


The  chamber  temperature,  Tc>  is 


(39) 


The  exit  temperature  (37)  is 


T 


C 


4 


The  exit  pressure  (36)  is 


P 

e 


C 


3 


(40) 


(41) 


where  Pc>  the  chamber  pressure,  is  usually  a  specified  engine  charac¬ 
teristic. 


The  engine  thrust  coefficient  as  evaluated  from  the  previously 
derived  parameters  is  (Hill  et  al  (1963)) 


jl±_L 


=  P-  (^-V  1 

>1  V-  1  \v+  1/ 


ft) 


and  the  maximum  thrust  coefficient  is 


max 


/  2  \  -  1 
F7'  v  +  1  / 


-  1 
V 


-T  -n^  (42) 


(43) 


At  this  point,  a  check  may  be  made  on  the  calculations  by  comparing  the 
thrust  coefficient  as  calculated  from  expression  (42)  with  that  derived 
from  specified  engine  characteristics 


T 


CF  =  'P  A* 


(44) 


or  by  comparing  the  total  engine  thrust  calculated  from 


T  =  P  A* 
F  c 


(4  5) 


with  the  specified  engine  thrust.  The  throat  area  A*  and  expansion  ratio 
A^/A*  are  usually  specified  quantities. 

The  plume  drag  of  a  rocket  with  two  phase  (gas  and  solid)  flow 
should  be  calculated  using  the  thrust  due  to  gas  alone  since  the  solid 
particles  do  not  take  part  in  the  expansion.  The  speed  of  sound  in  the 
gas  along  governs  the  expansion  of  the  plume  at  large  pressure  ratios. 
The  frozen  (gas  along)  Mach  number  is 

V 

M  =  e.L,  (46) 

g  -JyR  Te 

The  Mach  number  is  considerably  lov  sr  than  the  equilibrium 
(including  solids)  Mach  number.  The  thrust  of  the  rocket  engine  due 

to  the  gas  phase  is  calculated  by  determining  the  chamber  conditions 
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corresponding  to  the  gas  along  exit  conditions,  y  and  M  .  From  gas 

8 

tables,,  the  following  pseudo  engine  characteristics  are  obtained 


(47) 


(48) 


The  chamber  pressure  and  throat  area  of  the  pseudo  engine  are  deter¬ 
mined  since  the  exit  pressure  and  exit  area  are  known  quantities. 

The  pseudo  engine  thrust  coefficient  is 


2\yl 

v 

v 

Y  ! 

+  V 

\P 

c  hp 

The  maximum  thrust  coefficient  is 


max 


:(^r) 


1/y-i 


The  rocket  thrust  (gas  alone)  is 
T 


g 


=  (aA 

V  *  / p  \  c/p  V 


(49) 


(50) 


(51) 


For  comparison,  the  gas  thrust  can  also  be  determined  from 
T  =  T  (1  -  \)  (52) 

O 


where 


mass  flow  of  solids 
total  mass  flow 


(53) 
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free  stream  ratio  of  specific  heats 
Jq  a  constant  equal  to  .  8  5 
D  the  plume  drag 

free  stream  density 
vehicle  velocity 
The  maximum  plume  radius  is 


Cjj  Drag  coefficient  of  plume  shape 

P  free  stream  pressure 
oo  r 

M  free  stream  Mach  number 
00 

The  distance  to  the  Mach  disc  is 


limiting  velocity  of  exhaust  gas  flow 
VCx  exhaust  gas  velocity  at  exit  planes 
Mex  exit  plane  Mach  number 
T  thrust  due  to  gas  phase  flow 
Y  ratio  of  specific  heats  for  gas  phase  flow 

2?  MC 
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4.  AVERAGE  PROPERTIES  IN  THE  AIR -EXHAUST  LAYifR 


With  the  contact  surface  positioned,  it  is  possible  to  derive 
simplified  properties  of  the  double  shock  layer  (air-exhaust  gas  layer 
on  the  peiiphery  of  the  plume)  using  the  inviscid  model  outlined  by 
Alden  et  al .  By  treating  the  layer  as  if  it  were  a  mixture  of  two  gases, 
instead  of  being  separated  by  a  contact  surface,  it  is  possible  to  estimate 
the  overall  thickness  of  the  layer  and  its  average  properties.  The  thick¬ 
ness  of  the  layer  may  be  expressed  in  non-dimensional  form  and  is  a 
function  only  of  position  along  the  dividing  streamline  (i.  e.  ,  contact 
surface),  the  properties  of  the  exhaust  gases  and  the  ratio  of  the  free 
stream  velocity  to  the  jet  limiting  velocity.  In  functional  form,  6^^ 


6ML 


(55) 


A  new  parameter,  appears  in  the  determination  of  the  overall 

thickness  of  the  air-exhaust  gas  layer.  It  is  interesting  to  note  that  this 
parameter  influences  the  thickness  of  the  double  shock  layer  but  does 
not  influence  its  position. 


4.  1  Properties  of  Air  and  Exhaust  Gas  Layer 

The  analysis  may  be  extended  to  provide  the  properties  of  the 
individual  layers,  (i.  e.  the  properties  of  the  air  layer  between  the  con¬ 
tact  surface  and  the  outer  shock,  and  the  exhaust  gas  layer  on  the  in¬ 
side).  Equations  may  be  written  for  the  total  thickness,  the  mass  flow 
in  the  air  layer  and  exhaust  layers  and  the  total  momentum. 
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There  are  four  equations  with  four  unknowns 

T  T  VA»,  ml  d  VJET,ML 

3 JET, ML  ’  °AIR,  ML  *  VM  an  VM 

The  left  hand  sides  are  all  known.  The  densities  on  the  right  are  ob¬ 
tainable  from  the  pressure  and  enthalpy  by  using  an  averaged  equation 
of  state  for  each  layer. 


Pair,  ml 


(y1)air, 


AIR,  ML 
ML  hAIR, ML 


(62) 


PjET,  ML 


JET, ML 
JET,  ML  hJET,  ML 


(63) 


where 


V 


ML 


hAIR,  ML  + 


AIR,  ML  _  _o_ 


(64) 


.  JET,  ML  _  VM 

~njET,ML  +  2  "  2 


(h°)jET,;, 

Solving  Equation  (56)  through  (59),  ’8  ^oun(*  *° 

represented  by  a  solution  to  a  cubic' equation: 


M 


L) 


rv  \  ^  /v  ^  ^ 

VAI.R,ML  l  (  At-B+C  l  AIR,  ML 


D 


M 


) 


E+F+G+H 


r  AIR,  ML  I  i.  ,  0 

\  VM  I  D 
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With  ,Ml/^m)  determined, 

now  be  obtained. 


solutions  for  the  other  unknotvns  may 


This  technique  has  been  used  to  determine  the  characteristics 
of  the  overall  and  individual  layers  for  a  typical  plume  (Vj^t  =  32  5, 

h^exit  =  3.  0  and  =  1.  2  5).  The  ratio  of  specific  heats  for  overall 

layer  was  assumed  to  be  the  average  of  the  air  and  exhaust  gas  layers. 


The  variation  of  the  total  thickness  of  the  air-exhaust  gas  layer 
with  forward  speed  ratio  is  shown  in  Figure  14.  A  comparison  of  the 
total  layer  thickness  (Figure  14)  with  the  exhaust  layer  thickness 
(Figure  15)  for  a  V^/V^  riltio  near  unity  indicates  the  individual  layers 
are  nearly  the  same  thickness.  For  larger  /V^  ratios,  the  exhaust 
layer  becomes  much  thicker  than  the  air  layer.  The  average  pressure 
in  the  air -exhaust  gas  layer  is  noted  in  Figure  16.  The  mass  flov»  in 
the  overall  layer  and  the  mass  flow  in  the  exhaust  gas  layer  are  noted 
in  Figures  17  end  18  respectively.  Fifty  percent  of  the  mass  flow 
exiting  from  the  rocket  engine  is  entrained  in  the  exhaust  gas  leyer. 
when  the  non-  dimensional  longitudinal  distance,  X  equals  1.2.  The 
average  static  enthalpy  of  the  overall  layer  and  the  density  in  the  ex¬ 
haust  layer  are  shown  in  Figures  19  and  20.  The  magnitude  of  the 
density  in  the  exhaust  gas  layer  relative  to  the  exit  plane  density  may 
be  determined  from  Figure  20  if  the  dynamic  pressure  qQ,  the  engine 
chamber  pressure  P.  and  expansion  ratio  are  known.  If  frozen  fl  >w  is 
assumed  from  the  engine  exhaust  plane  through  the  exhaust  gas  layer, 
the  average  electron  density  in  the  exhaust  layer  is  also  determined. 

In  addition,  since  the  average  temperature  in  the  exhaust  gas  layer  is 
known,  the  layer  collision  frequency  may  be  determined. 
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Calculations  were  performed  for  the  case  where  the  specific 
heat  ratio  for  the  overall  layer  is  evaluated  using  the  relative  mass 
flow  in  the  individual  layers  instead  of  assuming  an  average  of  the 
exhaust  gas  and  air  layer  specific  heat  ratio,  i.  e.  ; 

_  CA'.R  CPA1R  *  CJET  CPJET 
'  CAIR  CVAIR  +  CJET  CVJET 

or 


The  variable  specific  heat  ratio  calculations  produced  results  sub¬ 
stantially  similar  to  the  constant  y  calculations.  The  variable  y  calcula¬ 
tions  predicted  slig''  Ty  smaller  magnitude  for  overall  and  individual 
layer  thickness  and  slightly  larger  values  for  the  density  in  the  exhaust 
layer.  The  mass  flow  in  the  overall  and  individual  layers  was  unaffected 
by  the  method  used  tc  calculate  y. 


5.  NEAR -FIELD  SHOCK  WAVES 
PRODUCED  BY  ROCKET  PLUMES 


.t'/i 


The  theory  describing  the  far -field  flow  field  about  supersonic, 
axisymetric  bodies  was  developed  by  G.  B.  Whitham  (1952).  This 
theory  can  be  used  to  predict  the  location  and  strength  of  shock  wave 
and  the  plume  drag.  Experimental  agreement  with  the  theory  has  been 
obtained  by  DuMond  (1946)  with  bullets,  Carlson  (1959)  in  wind  tunnel 
investigations  and  by  Moglieri,  (1959)  and  Mullen,  ,  i956)  in  flight  tests 
of  full  scale  airplanes. 

The  Whitham  theory  is  a  modification  of  linearized  supersonic 
flew  theory  which  introduces  shock  waves  by  replacing  the  Mach  lines 
of  the  linear  theory  with  accurate  expressions  for  the  characteristics. 
The  basic  distributing  function  of  the  body  shape  remains  unchanged. 

The  formulas  predicting  the  location  of  the  front  shock  wave, 
the  pressure  discontinuity,  and  the  plume  drag  are 


(66) 


(.67) 


F2(y),d  y 


where  j3 


polytropic  index 


(68) 


vlindrical  coordinates:  x  is  measured  along 
oody  axis  downstream  from  the  nose,  r  is  the 
radial  coordinate  measured  perpendicular  to  x. 
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M 


free  stream  Mach  number 


P  pressure 

D  drag 


AP  =  P  -  P 

00 

The  basic  function  F(t])  is  the  shape  integral  of  the  plume  defined 


F^-h  j- ^  <«> 

o 

2 

where  S  is  the  cross-sectional  area,  irR  .  The  value  of  the. upper  limit, 
y  ,  is  determined  from  the  condition  that  the  integral 


F(n)  d  t| 


(70) 


be  a  maximum  at  y  =  y  .  A  necessary  condition  for  a  maximum  at 

y  =  yG  is 


dl_ 

dy 


F(y)  =  o 


(71) 


The  Whitham  theory  is  inapplicable  to  an  elliptical  body  because 
of  the  blunt  nose.  However,  experimental  work,  Carlson  ( 1959),  shows 
th't  the  far-field  pressures  are  essentially  the  same  as  those  for  a  pointed 
parabolic  body  of  the  same  fineness  ratio.  Thus  instead  of  using  the  blunt 
elliptical  plume  shape  specified  as 
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the  plume  will  be  described  as  a  parabolic  body  having  the  same  rmax 
as  the  ellipsoidal  plume. 


Integrating  Eq.  (5.4),  the  shape  integral  is  obtained  as 


F(y )  = 


The  parameter  in  the  shock  wave  pressure  and  position  equa¬ 
tions  is  given  as 


F(y)  dy  =  0.  92  K 


where  the  upper  limited  is  determined  as  yQ  =  0.  692L. 

The  accuracy  of  this  approach  can  be  determined  by  evaluating 
the  plume  drag  using  Whitham' s  theory  (Eq.  68)  and  comparing  the 
results  with  the  known  drag  determined  by  Hill  and  Habert  (1963).  The 
drag  obtained  from  Eq.  (68)  is 


D  =  0. 4  5tt  p  V  2  K  2 
roo  oo  x 

The  drag  computed  by  the  Whitham  theory  agrees  to  within  3  percent  of 
the  drag  computed  for  the  ellipsoid  shape  using  the  drag  coefficients 
evaluated  by  Hill  and  Habert  (1963).  The  comparison  is  made  at  =  10 
on  a  representative  trajectory. 

A  wind  tunnel  test  at  Mach  number  2.01  is  available  (Carlson  1959) 
which  gives  direct  comparison  between  the  experiments  with  an  ellipsoid 
shape  and  the  theory  for  the  equivalent  parabolic  body.  The  comparison 
ic  shown  in  Figure  21. 
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as 


The  plume  shock  wavj  position  and  strength  can  now  be  described 


x  -  0r  -  0.  69  L  =  r1/4JT83k 


aP 

P 


The  shockwave  overpressure  hac  been  calculated  for  the  Saturn 
and  Atlas  second  stage  rockets  and  is  shown  in  Table  5.  1.  The  radial 
distance  from  the  flight  path  to  the  point  of  ten  percent  overpressure  is 
given,  neglecting  atmospheric  attenuation  of  the  shockwave.  The  above 
expression  describes  the  shockwave  to  the  point  it  decays  into  an  acoustic 
wave.  Thereafter,  the  acoustic  disturbance  must  be  treated  by  taking 
into  account  the  variation  in  atmospheric  transport  properties  using  ray 
tracing  methods  or  specifying  the  technique  developed  by  Meyer  (1962). 
The  near-field  model  describing  shock  waves  of  strength,  aP/P^  >  0.  1, 
is  valid  for  distance  up  to  50-  100  km  from  the  plume  (see  Table  5.  1). 
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6.  SUMMARY  AND  CONCLUSIONS 


Theoretical  models  ere  d' scribed  which  provide  an  improved 
analytic  description  of  the  gas  dynamics  of  high  altitude  rocket  plumes. 
These  methods  are: 

1.  A  simple  analytic  description  of  the  exhaust  flow 
from  rocket  motors  into  the  region  bounded  by  the 
inner  shock  wave  and  Mach  disc. 

2.  A  model  for  the  two  phase  flow  in  the  exhaust  of 
solid  propellant  rocket  vehicles  and  the  effect  of 
such  flow  on  the  resulting  plume  characteristics. 

3.  An  analytic  description  of  the  flow  properties  in 
the  air-exhaust  gas  layer  on  the  periphery  of  the 
rocket  plume. 

The  results  described  in  this  report  in  conjunction  with  those 
previously  developed  by  Hill  et  al  and  Alden  et  al  allow  the  overall 
features  of  plumes  from  both  liquid  and  solid  propellant  rocket  vehicles 
to  be  described. 

In  addition  the  structure  of  plumes  has  been  described  in  suffi¬ 
cient  detail  to  allow  an  ionization  model  for  high  altitude  rocket 
plumes  to  be  constructed.  Such  a  model  is  described  in  an  Addendum 
to  this  report. 
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Figure  1.  Structure  of  high  altitude  plume 
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Figure  2.  Schematic  flow  pattern  of  nozzle  exhausting 
into  a  vacuum. 
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Figure  8.  "Exact"  and  approximate  calculations  of  density 
decay  along  plume  axis. 
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Figure  11.  Coordinates  for  calculation  of  idealized 
plume  surface. 
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PLUME  SCALING  PARAMETER 
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iguxe  21.  Comparison  of  Whitham  theory  with  experiment 
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